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1.  INTRODUCTION 


Tlie  Combined  Release  and  Radiation  Effects  Satellite  (CRRES)  is  a  spin  stabilized  spacecraft  in  a  near 
equatorial  eccentric  orbit,  with  perigee  around  350  kilometers  and  apogee  near  33,500  kilometers.  The 
CRRLS  attitude  is  essential  for  conversion  of  magnetometer  data  to  geophysical  coordinates  and  for  the 
interpretation  of  solar  panel,  photometer  and  other  experiments  on  board,  llie  satellite  is  equipped  with 
a  sun  sensor,  a  horizon  sensor  and  a  low  resolution  magnetometer  for  use  in  attitude  determination. 
With  proper  treatment,  these  instruments  are  capable  at  present  of  resolving  the  attitude  of  the  spacecraft 
to  better  than  2®. 

The  CRRES  Attitude  Determination  Program  (CADP)  was  developed  by  Space  Applications 
Corporation  [1989]  and  delivered  in  July  1989,  almost  a  year  before  launch  of  the  satellite.  Ihe  attitude 
determination  and  modeling  routines  were  designed  for  nominal  operations  of  the  s|).itecraft.  Alllioiigh 
the  actual  definition  of  ‘nominal  operations’  has,  to  our  knowledge,  never  been  clearly  set  down,  one  can 
infer  the  requirements  from  the  design  specifications  of  the  original  CADP.  These  limitations  include 
(a)  pu’e  rotation  around  the  spacecraft  z-axis,  (b)  adequacy  of  the  averaged  spin  rate  calculated  from  sun 
sensor  data  for  linear  representation  of  the  phase  during  each  period  or  segment  delineated  for  modeling 
attitude,  and  (c)  high  accuracy  in  the  engineering  magnetometer. 

Upon  delivery,  Radex  undertook  the  responsibility  for  integration  of  the  system  into  the  overall  Orbital 
Data  Processing  (ODP)  program  and  for  monitoring  the  performance.  Before  launch,  it  became 
apparent  from  simulated  data  that  the  linear  model  of  the  phase  would  be  severely  limiting  for  some 
applications.  Even  with  a  constant  spin  rate,  propagation  of  errors  in  spin  axis  position  and  natural 
motior  of  the  sun  in  the  defining  coordinate  system  gave  rise  to  a  drift  in  the  calculated  spin  rate  and 
required  two-hour  segmentation  with  approximately  one-half  degree  discontinuities  at  segment 
boundaries.  To  remove  this  limitation,  an  attitude  refitting  program  was  developed  [McNeil,  1990]  to 
use  the  calculated  phase  to  find  a  spin  rate  correction  term.  Originally,  this  was  intended  to  be  used  only 
to  reprocess  the  attitude  for  use  in  magnetic  field  despinning.  However,  shortly  after  launch,  it  became 
apparent  that  specifying  a  constant  spin  rate  even  under  normal  operation  would  lead  to  unacceptable 
errors.  It  was  also  clear  that  the  calculation  of  attitude  from  magnetometer  data  was  inadequate  except 
for  short  periods  around  perigee. 

In  response  to  these  problems,  the  attitude  refitting  system  was  enhanced  to  include  a  linear  and 
eventually  a  quadratic  spin  rate  during  sunlit  periods  and  a  quadratic  spin  rate  during  eclipse.  Also,  a 
global  spin  axis  pointing  direction  was  determined  using  sun  and  horizon  sensor  data  alone,  lliis 
software  has  been  used  exclusively  in  the  generation  of  the  CRRES  attitude  model.  Jlie  modeling  is 
done  interactively  from  intermediate  results  of  the  original  CADP,  and  the  model  coefficients  generated 
by  the  CADP  are  discarded.  Integration  of  the  revised  model  into  the  batch  processing  of  attitude  has 
never  been  attempted  for  a  variety  of  reasons.  First,  times  for  attitude  adjustments,  spin  rate  changes 
and  other  anomalies  are  not  generally  available  at  run  time  and  determination  of  these  events  is  quite 
easily  accomplished  from  inspection  of  the  intermediate  attitude  data.  Also,  telemetry  anomalies  and 
missing  data  can  be  more  easily  dealt  with  by  interactive  segmentation.  Finally,  the  interactive  fitting 
requires  little  extra  time  for  all  .  ut  the  most  difficult  orbits,  and  the  system’s  presentation  graphics 
provides  a  high  degree  of  quality  control. 
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'Fhis  document  is  meant  to  serve  several  purposes.  At  present,  the  only  formal  documentation  of  the 
attitude  determination  algorithms  is  a  design  specification  written  in  Program  Design  Language 
(commonly  called  PDL).  In  this  form,  it  is  extremely  difficult  to  follow  the  calculations  from  beginning 
to  end.  Thus  it  seems  valuable  for  future  reference  to  describe  the  calculations  in  more  concise 
mnthematical  terms,  albeit  that  most  are  relatively  standard.  Second,  there  is  a  need  to  document  the 
changes  made  to  the  system  since  delivery  and  especially  the  attitude  modeling  software,  lliird,  it  seems 
valuable  to  present  some  firm  accuracy  estimates  for  the  attitude.  Finally,  this  seems  a  good  medium  for 
the  presentation  of  some  of  the  idiosyncracies  of  the  CURES  attitude  behavior  to  serve  as  a  guide  in 
thinking  about  attitude  software  development  for  similar  missions. 

In  Section  2,  the  attitude  instrumentation  for  CRRES  is  briefly  described  and  the  algorithms  for 
calculation  of  spacecraft  spin  axis  pointing  direction  and  spin  phase  from  simultaneous  observations  are 
presented.  Next,  we  discuss  the  specific  problems  encountered  in  the  generation  of  the  attitude  model 
and  the  solutions  presently  in  place.  Finally,  we  describe  the  anomalies  one  might  encounter  in  the  use 
of  the  CRRES  attitude  model  and  discuss  the  magnitude  of  uncertainty  to  be  expected. 


2.  ATTITUDE  DETERMINATION 


lliis  section  describes  the  instruments  and  algorithms  used  to  determine  an  ‘instantaneous’  attitude  point 
for  the  satellite.  The  CRRES  attitude  is  defined  in  Earth  Centered  Inertial  (ECI)  coordinates.  In  ECI, 
the  z-axis  is  collinear  with  the  Earth’s  rotation  axis  and  the  x-axis  points  toward  the  vernal  equinox  in 
the  ecliptic  plane.  Tlie  attitude  is  defined  by  the  right  ascension  and  declination  of  the  spin  axis  (the 
CRRES  z-axis)  and  by  a  phase  which  is  measured  in  the  spacecraft  spin  plane  relative  to  passage  of  the 
spacecraft  x-axis  through  ihe  ECI  xy-plane  traveling  northward.  The  attitude  ‘point’  is  comprised  of  the 
rigiit  ascension  a  and  declination  5  of  the  spin  axis,  an  instantaneous  spin  rate  co,  and  a  time  t^  at  which 
the  spin  phase  was  zero. 


2.1  INSTRUMENTA'IION 


In  this  discussion,  the  attitude  instrumentation  aboard  the  satellite  will  be  described  only  to  the  extent 
needed  to  understand  the  data  obtained  from  them.  Details  of  the  instruments  themselves  have  been 
given  elsewhere  [Ball,  ]986a,b,c].  'Hie  CRRES  attitude  instruments  are  shown  schematically  in  Figure 
1.  The  sun  sensor  consists  of  a  fan  located  parallel  to  the  spin  axis  on  the  negative  x  side  of  the 
spacecraft.  Each  time  the  sun  passes  through  the  fan,  the  angle  between  the  spin  axis  and  the  sun  is 
determined  to  an  accuracy  of  +A  0.5°.  Also,  the  time  of  the  crossing  of  the  sun  through  the  spacecraft 
xy-plane,  tggi  >s  determined.  'Lhis  time  is  inserted  into  telemetry  in  a  rather  complex  way,  as  are  the  other 
event  times  discussed  below.  The  details  of  this,  which  will  be  avoided  here,  are  ava.Iable  in  Bah 
Aerospace  documentation  (1985]. 

The  attitude  magnetometer,  referred  to  as  the  engineering  magnetometer  in  what  follows  to  differ .jntiate 
it  from  the  science  magnetometer,  is  tri-axial  and  aligned  with  the  spacecraft  axes  to  better  than  0.2°. 
It  samples  the  field  approximately  once  per  second  in  two  gain  states,  one  with  resolution  of 
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approximately  4S0  nT  and  another  with  approximately  74  nT  resolution.  Also,  the  time  t^^g  at  which 
the  field  along  the  magnetometer  x-axis  is  zero  is  determined  and  inserted  into  telemetry. 

The  Earth  horizon  sensor  provides  the  times  that  the  50%  radiance  point  of  the  Earth’s  IR  horizon 
crosses  the  sensor’s  field  of  view.  These  are  separated  into  leading  edge,  tj^LE'  trailing  edge,  tj^5.f£. 
Contrary  to  the  depiction  in  Figure  1,  the  horizon  sensor  is  now  permanently  fixed  in  the  spacecraft  spin 
plane  and  is  not  steerable. 


2.2  TELEMETRY  PROCESSING 


lire  first  step  in  attitude  determination  for  CRRES  is  the  unpacking,  time-tagging  and  conversion  to 
physical  units  of  the  raw  attitude  telemetry'  values.  The  data  file  used  is  a  reduced  version  of  the  total 
telemetry  containing  data  from  the  three  attitude  instruments  only.  Each  Master  Frame  (MF)  in  this 
file  has  been  time  tagged  with  the  Universal  "lime  (UT)  at  the  start  of  the  master  frame  by  the  Master 
Frame  Formatting  (MFF)  program.  'ITie  processing  of  the  sun  sensor  data  requires  the  conversion  of 
the  measured  sun  declination  angle  from  its  ‘grey  scale’  code  to  degrees  [McNeil  and  Mclnerney,  1988] 
and  the  calculation  of  the  ‘absolute’  UT  of  the  sun  hit. 

Time  tagging  of  the  sun  pulse  is  identical  to  that  of  the  horizon  sensor  and  magnetometer  zero  aossing 
data.  Tliese  are  derived  from  a  local  time  counter  with  a  resolution  of  64  kHz.  The  local  clock  is  reset 
to  zero  at  the  master  frame  sync  pulse.  The  pulse  occurs  1.5  ms  before  the  start  of  each  master  firame. 
When  a  sensor  event  occurs,  the  value  of  the  local  clock  is  latched  by  the  appropriate  holding  register. 
This  value  remains  constant  in  the  register  until  the  next  such  event  occurs  and  is  dumped  twice  per 
major  frame.  In  order  to  determine  the  actual  time  of  a  sensor  event,  one  first  examines  the  values  of 
successive  data  to  see  when  a  change  in  the  value  has  taken  place.  One  must  next  determine  if  the  local 
clock  time  refers  to  the  present  or  the  previous  master  frame.  For  example,  the  first  of  the  two  sun 
pulses  is  read  out  in  word  62  of  subframe  1.  'Iliere  are  32  subframes  per  master  frame  (numbered  0 
through  31)  and  256  words  per  subframe  (numbered  0  through  255).  The  readout  time  of  the  first  sun 
pulse,  relative  to  the  MF  time,  is 


t,  =  0.128  (1  +  62/256)  sec. 

at  the  nominal  rate  of  4.096  seconds/MF.  'Hie  value  of  the  local  time  clock  is  found  from 

t^  =  .0000625N  -  .0015  sec. 

where  N  is  the  clock  count,  and  the  difference  between  MF  sync  pulse  and  the  actual  start  of  the  MF 
has  been  subtracted  to  make  t^  relative  to  the  start  of  the  MF.  If  t^  is  less  than  t,  then  the  readout  must 
have  pertained  to  the  current  master  frame,  and  the  time  of  the  crossing  is  given  by 

*SS  “  *MH  *e 

where  t55  is  the  assigned  time  of  the  current  master  frame.  Otherwise,  the  time  is  given  by 

*ss  ~  *MF  *c '  ■*•096  sec. 
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Additionally,  the  horizon  sensor  has  an  additional  built-in  delay  time  tj  which  must  be  added  to  tlie  time 
as  calculated  above.  'Fhis  time  is  taken  to  be  645.55  ms  for  2  rpm  mode  and  51.54  ms  for  10  rpm  mode. 

Finally,  the  analog  readings  of  the  engineering  magnetometer  are  converted  to  the  magnetic  field  values 
along  the  three  spacecraft  principal  axes. 


2.3  OBSERVATION  SETS 


From  the  calibrated  telemetry  data,  a  series  of  observation  sets  is  next  formed.  Each  of  these  constitute 
data  taken  during  one  spin  period.  Some  of  these  sets  contain  pairs  of  observations  from  which  the 
attitude  Gin  be  cilculated  in  several  different  ways.  Others  contain  insufficient  data  to  Gilculale  it  even 
once.  Observation  sets  are  defined  from  one  sun  hit  to  the  next  or,  in  eclipse,  from  one  crossing  of  the 
magnetic  field  zero  to  another.  For  the  CRRES  instruments,  the  possible  combination  of  measurements 
is  given  in  Table  1. 


Table  1.  Possible  CRRES  Observation  Sets 


# 

Sensors 

When  Applicible 

I 

SS  Alone 

Never  (0%) 

2 

HS  Alone 

Never  (0%) 

3 

SS,HS 

Never  (0%) 

4 

Mag  Alone 

Typical  Eclipse  (8%) 

5 

SS,  Mag 

Full  Sun,  Earth  not  in  HS  FOV  (75%) 

6 

HS,  Mag 

Occasionally  during  Eclipse  (2%) 

7 

SS,  HS,  Mag 

Full  Sun,  Earth  in  HS  FOV  (15%) 

The  first  four  of  these  events  are,  as  indicated,  not  applicable  to  the  way  in  which  the  CRRES  data  is 
currently  processed.  Had  the  use  of  the  magnetic  field  been  optionally  limited  to  the  regions  very  close 
to  the  earth,  the  majority  of  the  data  would  have  been  of  type  1.  For  reasons  treated  in  more  detail  later, 
we  have  chosen  not  to  limit  the  use  of  the  magnetometer  at  this  level,  'llie  percentages  given  above 
apply  roughly  to  typical  orbits  during  the  first  six  months  or  so  of  the  mission.  As  can  be  seen,  only  a 
relatively  small  portion  of  the  data  contains  earth  sightings. 


2.4  SPIN  RATE 


From  comparing  the  time  of  each  observation  set  \rith  that  of  the  next,  assuming  both  are  defined  by 
a  sun  or  magnetometer  zero  cross  pair,  the  instantaneous  spin  rats  is  determined.  The  calculation  for 
sun  crossing  pairs  is  straightforward.  That  for  magnetometer  zero  cross  pairs  is  equally  so  except  for 
a  correction  term  added  later  that  deserves  mention.  The  correction  can  be  carried  out  only  when  an 
estimate  is  available  for  the  spin  axis  pointing  direction,  z^.  In  the  CRRES  package,  the  estimate  used 
is  that  of  the  first  guess  attitude,  the  calculation  of  which  is  explained  below.  The  correction  goes  as 
follows.  The  ECI  unit  vectors  of  the  model  field  at  the  first  zero  cross,  bj  and  the  second  zero  cross, 
b2  are  found.  Taking  the  cross  product  of  each  with  the  spin  axis  gives  two  vectors,  Cj  and  Cj  which  are 
equivalent  in  phase  to  the  projection  of  the  magnetic  field  into  the  spin  plane  at  the  times  of  the  zero 
crossings.  Taking  the  phase  angle  betw'een  Cj  and  C2  gives  a  correction  term  for  the  spin  rate  deduced 
from  the  two,  representing  the  motion  of  the  background  field  betw'een  the  two  zero  crossings. 
Specifically,  if  is  the  spin  rate  found  from  subtracting  times  for  two  successive  zero  crossings,  then 
the  corrected  spin  rate 


<>»nug  = 

where  6  is  the  rotation  angle  from  Cj  to  €2-  This  rotation  angle  is  calculated  from  the  law  of  cosines. 
Since  this  particular  transformation  wall  come  up  again  in  what  follows,  we  define  the  rotation  angle 
about  the  vector  a  from  the  projection  of  vector  fa  into  the  plane  perpendicular  to  a  to  the  projection 
of  the  vector  c  into  the  same  plane,  A(a,b,c),  on  the  inten-al  [0,2jti  as  follows. 


cost|r  = 


6'c  —  (a*c)(6‘c) 


(1) 


A(a,b,c)  =  cos**(costJf)  c-(axb)  >  0 

A(a,b,c)  =  2ic  -  cos‘*(cost|»)  c-(axb)  <  0 

'Ilie  complement  operation  above  simply  puts  the  angle  in  the  proper  quadrant.  In  the  calculation 
above,  0  =  A  (2cji,Cj,C2).  Whenever  a  proper  spin  rate  cannot  be  calculated  for  an  observation  set,  the 
so  called  telemetered  spin  rate  is  substituted,  'llic  telemetered  spin  rate  is  calculated  on  board  from  an 
average  of  the  last  sixty  sun  hits. 


6 


2.5  FIRST  GUESS  A'lTITUDE 


For  each  pair  of  observations  in  each  set,  two  distinct  spin  axis  pointing  directions  can  be  calculated. 
In  order  to  determine  which  one  is  correct,  one  must  do  some  sort  of  survey  of  a  large  number  of  points 
taken  from  different  positions  in  an  orbit  or  from  different  sensors.  The  false  solution  will  drift  as  the 
objects  sensed  change  orientation  in  the  spacecraft  frame  while  the  true  solution  will  remain  relatively 
constant.  For  CRRES,  this  is  done  by  calculating  the  spin  axis  pointing  direction  z  from  the  first  one- 
hundred  observations  using  the  sun  sensor  data  plus  tj^gm,  JnsTE  ^mag*  calculated  as 

follows. 


Sun  Sensor /Magnetometer  Attitude.  Using  the  times  of  sun  hit  and  magnetic  field  zero  cross, 
we  find  the  rotation  angle  between  the  sun  and  the  spacecraft  x-axis  at  the  time  of  the  zero  cross 

t  =  2n«(t„,g  -  tss) 

The  angular  separation  between  the  x-axis  and  the  sun  at  the  zero  cross  is  given  by  the  law  of  cosines, 
which  reduces  here  to 


costi  =  sinPcos<t 

where  P  is  the  sun  aspect  angle.  We  find  frvo  possible  ECI  pointing  directions  for  the  spacecraft  x-axis 
using  the  ECI  model  reference  vectors  for  the  sun  s  and  magnetic  field  b  by  taking  the  interception  of 
frvo  cones  around  s  and  b  with  half-vertex  angles  t)  and  n/2  respectively. 

A  general  method  f  Wertz,  1986]  is  as  follows.  Suppose  the  vectors  e  and  s  have  aspect  angles  and 
P  respectively.  This  is  shown  schematiailly  in  Figure  2.  We  lake 


cos  p-«*«  costi 
1  - 


COST] -e‘s  cos  P 


(2) 


z=± 


l-zcosP-y  cost! 


l-(e'8f 


'Ilien,  forming 


c  =  s  X  e 
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Figure  2.  Determination  of  degenerate  attitude  vector  using  two  known  vectors  and  two 
arc  length  measurements.  From  Wertz  [1986]. 


we  find  the  interception  vectors  from 


a  =  as  +  j/e  +  2C 

We  refer  to  this  calculation  in  what  follows  as  the  function  S(e,ti,s,P).  In  the  calculation  of  the  auxiliary 
reference  vectors  above,  we  have 


hj  2  =  S(s,Ti,b,7t/2) 

From  here,  two  possible  spin  axis  pointing  directions  are  calculated  by  finding  the  interception  of  the 
sun  vector  s  with  half-vertex  angle  P  and  the  vectors  hj  2  again  with  vertex  angle  it/2.  'I  his  gives  a  total 
of  four  calculated  spin  axis  pointing  directions  in  all. 

2ia,b  =  3(hi,it/2,s,p) 

Z2a,b  =  S(h2,K/2,s,p) 

From  each  pair,  say  Zjg  j,,  there  will  be  one  that  does  not  reproduce  the  measured  rotation  angle.  By 
examining  the  rotation  angle 


a  =  A(zj3,s,h,) 

in  comparison  to  the  measured  angle  $,  we  can  tell  which  one  to  accept,  thus  reducing  the  total  number 
of  solutions  to  two. 

Sun  Sensor/Horizon  Sensor  Attitude.  For  sun  with  earth  leading  or  trailing  edge  calculations, 
a  spherical  earth  is  assumed  for  first  guess  attitude.  The  apparent  radius  of  the  earth’s  CO2  horizon  is 
calculated  from 


p  =  sin'^[(r^  +  rco2)/rsat] 

where  r^  is  the  mean  radius  of  the  earth,  r(-Q2  is  the  altitude  of  the  CO2  horizon,  taken  as  40  km,  and 
r^gj  is  the  length  of  the  satellite  to  earth  vector,  llie  angular  separation  between  the  reference  vector 
(horizon  sensor)  and  the  sun  is  again  found  from  the  law  of  cosines, 

COST)  =  cospcosy  +  sinPsinycos^* 

where  p  is  the  sun  aspect  angle,  y  is  the  horizon  sensor  cant  angle  (which  is  it/2  in  practice  but  variable 
in  the  code)  and  «!>  is  the  rotation  angle  between  sun  and  earth  hits.  We  find  two  possible  horizon 
sensor  vectors  from  the  ECI  sun  vector  s  and  the  satellite-to-earth  vector  e  by  intersecting  the  cones  ti 
away  firom  s  and  p  away  from  e. 


hj  2  =  S(s,Ti,e,p) 

ITien,  the  four  possible  spin  axis  positions  are  found  from 

2ta,b  =  2(s,p,hi,y) 
22a,b  =  S(s,p,h2,y) 
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and,  as  before,  one  of  each  pair  is  eliminated  by  examination  of  the  rotation  angle. 

First  Guess  Histogram.  From  the  Orst  100  attitude  points  calculated,  a  histogram  is  constructed 
to  determine  the  first  guess  attitude.  Originally,  more  than  one  first  guess  interval  was  possible  during 
a  single  orbit.  Additional  first  guess  intervals  were  added  at  times  of  attitude  adjusts.  However,  lack  of 
knowledge  of  mission  events  prior  to  execution  of  the  attitude  calculation  forced  this  approach  to  be 
abandoned.  This  could  have  been  a  serious  problem  except  for  the  fact  that  the  pre-adjust  and  post¬ 
adjust  attitudes  are  close  enough  so  that  the  same  first  guess  works  equally  well  for  both.  This  is  a  result 
of  the  CRRES  attitude  itself.  Since  the  spin  axis  points  nearly  at  the  sun  and  since  the  sun  drifts 
counterclockwise  looking  down  the  ECI  z-axis,  the  spacecraft  spin  axis  will  ‘trail’  the  sun.  The  sun  drifts 
away  from  the  spin  axis  to  the  maximum  allowed  angle  of  around  15°.  Then  a  series  of  precession 
maneuvers  are  performed  to  push  the  spin  axis  back  to  about  5°.  Since  the  two  possible  solutions  using 
sun  data  are  on  opposite  sides  of  the  sun  and  since  the  post-adjust  spin  axis  attitude  is  on  the  same  side 
as  the  pre-adjust  attitude,  the  post-adjust  true  solution  is  closer  to  the  first  guess  for  the  pre-adjust 
attitude  than  is  the  false  solution.  It  is  therefore  selected  using  the  pre-adjust  first  guess. 

The  solutions  are  binned  by  an  angle  defined  by  the  rotation  angle  between  the  attitude  solution  z 
around  the  mean  sun  vector  for  the  orbit  s  fi-om  a  vector  close  to  the  ECI  z-axis.  Because  the  CRRES 
spin  axis  points  nearly  toward  the  sun,  this  scheme  works  nicely  for  attitude  points  determined  from 
horizon  sensor  data.  However,  because  the  magnetic  field  is  nearly  perpendicular  to  the  satellite-to-sun 
vector,  true  and  false  solutions  derived  from  magnetic  field  zero  crossings  alone  tend  to  form  equal 
histograms  on  both  sides  of  the  sun.  Originally,  the  program  took  the  first  100  points  only,  not  generally 
including  horizon  sensor  data.  Because  of  this,  ambiguous  histograms  and  incorrect  first  guesses  often 
resulted.  This  situation  is  quite  serious  since  it  could  result  in  incorrect  choice  of  the  spin  axis  pointing 
direction  throughout  the  orbit.  We  found  a  solution  by  restricting  the  first  guess  histogram  calculation 
to  observation  sets  containing  earth  limb  data.  A  comparison  of  histograms  with  and  without  this 
restriction  is  given  in  Figure  3  and  Figure  4. 

First  guess  determination  using  magnetic  field  data  alone  is  not  able  to  reliably  differentiate  between  the 
true  and  false  solutions.  In  this  case,  the  most  populated  bin  would  indeed  have  been  the  correct  one, 
but  in  other  cases  this  is  not  true.  Even  this  aise  is  rather  too  close  for  comfort  with  the  most  populous 
bin  only  about  20%  more  frequent  that  one  far  from  the  actual  solution,  lliese  histogram  plots  are 
produced  from  the  attitude  modeling  and  quality  control  program  developed  by  Radex,  Inc.  Plots  like 
this  are  viewed  for  each  orbit  processed  in  order  to  assure  an  adequate  first  guess. 


2.6  A'lTI  rUDE  S'l’A'I  E  VECI  ORS 


Once  the  first  guess  for  the  spin  axis  pointing  direction  has  been  determined,  each  pair  of  useful 
observations  in  each  of  the  observation  sets  is  used  to  calculate  a  two  fold  degenerate  ‘instantaneous’  spin 
axis  pointing  vector  at  the  time  of  the  observation  set.  'Ilie  first  guess  is  used  to  choose  the  correct 
attitude  from  this  set  and  the  attitude  state  vector  for  the  observation  set  is  calculated  from  a  weighted 
average  of  all  calculations  made  within  that  set.  ITie  weighting  factors  are  three  in  number  and  are 
defined  in  the  CADP  input.  'Fhere  is  one  for  the  sun  sensor,  one  for  the  magnetometer  and  a  third  for 
the  horizon  sensor.  The  average  state  vector  is  calculated  as  follows. 
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FIRST  GUESS  HISTOGRAM  --  ORPTIT  0590 
INTERVAL  1  OF  1 

CHOSEN  BIN  49.  R. A.  8.4  DEC  -2.4 


Figure  3.  First  guess  histogram  from  the  original  CRRES  Attitude  Determination 
Program. 
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NUMBER  OF  OCCURENCES 


FIRST  GUESS  HISTOGRAM  —  ORBIT  0590 
INTERVAL  1  OF  I 

CHOSEN  BIN  49.  R.A.  8.4  DEC  -2.5 


Figure  4.  First  guess  histogram  obtained  by  restriction  of  the  observation  sets  included 
to  those  containing  earth  limb  data. 
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NUMBER  OF  OCCURENCES 


(3) 


^ave "  E  ^SS/HS  “^SS'^ag  S  ^SSMag 

(SS/liS  pairs)  (SS/Mag  pairs) 

SS/Mag  Pairs.  These  are  calculated  identically  to  those  calculated  in  Section  E.  'ITie  solution 
nearest  the  first  guess  is  saved  for  mixing  with  the  other  solutions. 

SS/HS  Pairs.  There  are  two  spin  axis  calculations  carried  out  with  these  data  iur  observation 
sets  with  sun  seirsor  and  horizon  sensor  data.  They  are  performed  in  a  manner  similar  to  that  described 
in  Section  '£  with  the  exception  that  an  iterative  correction  is  made  for  the  earth’s  oblaleness.  First,  the 
algorithm  in  Section  E  is  performed  with  the  spherical  earth  estimate  for  the  horizon  sensor  to  nadir 
angle  p.  'Ilien,  using  the  ailculated  value  of  h,  or  hj,  the  choice  being  deduced  from  selection  of  the 
proper  Zj  or  Zj,  a  revised  value  of  p  is  calculated.  'Ihe  algorithm  for  finding  this  value  is  given  in  Wertz 
(19^)  pp.98-101  and  will  not  be  repeated  here.  The  process  is  repeated  a  few  times  to  insure 
convergence. 

IJS/Mag  Pairs.  Although  there  is  a  rather  elaborate  calculation  of  attitude  from  horizon  sensor 
leading  and  trailing  edge  plus  magnetometer  zero  crossing  included  in  the  system,  this  is  used  only  when 
neither  of  the  two  combinations  above  can  be  found.  The  details  of  this  calculation  can  be  found  in 
Wertz  pp.366-370.  This  means  that  the  calculation  is  used  only  during  eclipse  when  earth  horizon  hits 
occur.  Tliis  happens  very  sporadically  and,  as  such,  is  not  pait.cularly  useful  for  routine  processing. 
Aside  from  this,  it  has  been  our  experience  that  this  computation  almost  nev^r  succeeds  in  producing 
reasonable  values.  This  is  most  likely  because,  as  we  will  see,  the  magnetic  field  data  is  not  veiy  useful 
beyond  a  few  thousand  kilometers  altitude. 

Spin  Phase  Zero  Time.  When  the  possible  spin  axis  attitude  vectors,  enumerated  above,  have 
been  calculated  and  mixed  according  to  the  weights,  the  time  of  the  zero  phase  crossing  for  each 
observation  set  is  calculated.  At  the  time  of  zero  phase,  the  spacecraft  x-axis  lies  in  the  ECI  x-y  plane 
and  is  rising,  by  definition.  Tliis  vector  is  given  by 

h  =  k  xz 

where  k  is  the  ECI  unit  vector  and  z  is  the  calculated  spin  axis  vector  for  the  observation  set.  Tire 
rotation  angle  from  the  sun  vector  s  to  the  x-axis  vector  h  is  calculated  by 

=  A(z,s,h) 

since  we  know  the  time  of  the  sun  crossing,  tgg  we  can  find  the  time  of  the  zero  phase  crossing  from 

tjj  “  tg<;  4" 

where  co  is  the  spin  rate,  llrere  is  an  equivalent  a.lculation  involving  the  magnetometer  zero  cross  time 
that  is  used  when  the  sun  is  not  visible.  Here,  the  model  field  vector  b  is  substituted  for  the  sun  vector 
s  and  the  magnetic  field  crossing  time  for  the  sun  crossing  time  tgj. 
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2.7  CADP  SEGMENTATION  AND  FITTING 


Although  the  attitude  model  produced  by  CADP  is  discarded  in  the  present  processing  scheme,  it  is 
valuable  to  examine  the  details  of  this  model  in  order  to '  ■’.s  why  the  results  were  not  acceptable.  With 
a  set  of  state  vectors  for  an  orbit  in  hand,  the  CADP  d; .  x  ",  orbit  into  segments.  A  segment  occurs 
at  the  beginning  and  end  of  eclipse,  at  each  mission  it  to  the  program  and  in  lieu  of  these,  at 

equal  intervals  of  predefined  duration,  set  cunently  ir  ■  For  each  segment,  the  state  vector 

times  are  normalized  to  the  interval  [-1,1]  via  the  u>r'<:jl!* 


t-t, 

X  =  2 - 1 


<2- 


'i 


(4) 


'Hiey  are  then  fit  by  standard  IMSL  routines  [IM  'L,  1987]  tc  Chebyshev  polynomials.  The  first  few  of 
these  polynomials  are 

go(X)  =  1  g,(X)  =  X 

g2(X)  =  2X2-1  e3(X)-4X3-3X 

g4(X)  =  8X^-  8X2 

'0)6  fitting  of  the  right  ascension  a  and  declination  8  is  carried  out  by  increasing  ,he  order  of  the  fit 
until  the  change  in  the  residual  by  increasing  the  order  is  ao  .more  than  a  preset  parar.ie'.er  e.  The  spin 
rate  is  fit  to  the  same  polynomial  expansion,  however,  the  order  is  never  allowed  to  go  beyond  a 
constant.  This  is  m  .essaiy  because  the  spin  rate  is  the  derivative  of  the  spin  phase,  which  is  actually  t'.ie 
attitude  parameter  of  importance.  Thi'i  ,i!iase  is  evaluated  by  the  Attitude  Agency  Module  (AGMOD) 
from 


(j)(t)  =  «j)o  +  co(X)(t-ti) 

where 

=  "ogoW  +  +  "zSiW 

Were  the  spin  rate  to  change  within  a  segment,  and  were  the  change  to  be  modeled  by  even  a  linear  term 
in  &>(X),  this  would  be  incorrect  since  the  phase  is  actually  given  by 

t 

+  J  (5) 


We  will  return  to  this  complication,  and  its  resolution,  in  the  section  on  attitude  modeling. 
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In  any  case,  with  a  constant  spin,  rate,  the  phase  at  the  beginning  of  each  segment  is  ailculated  in  the 
CADP  as  follows.  For  each  pha^e  zero  time,  the  initial  phase  is  calculated  by  extrapolating  backwards 
to  the  st.'i.rt  of  the  segment. 


-  ‘i) 

'.liese  are  averaged  to  give  the  value  used  for  the  ',;hase  at  the  start  of  the  segment,  (j)Q. 


3.  CRRES  ATTITUDE  CHARACTERISTICS 


In  order  to  understand  the  problems  encountcicd  in  modeling  the  CRRES  attitude  an  :  .he  reasons  why 
the  batch  processed  model  was  abandoned,  it  is  necessary  to  point  out  some  of  the  <  ..tinctis'e  features 
of  the  CRRES  attitude.  The  CRRES  mission  j  divided  into  two  separate  classes  oS  experiments,  the 
Geosynchronous  Transfer  Orbit  (GTO)  experiments  and  the  Low  Altitude  Satellite  St.idies  of 
Ionospheric  Irregularities  (LASSII)  experiments.  Tire  GTO  phases  of  the  mission  constitute  >95%  of 
the  data  gathering  periods.  The  LASSII  periods  are  almost  always  less  than  15  minutes  in  duration  and 
present  less  of  a  problem  for  attitude  determination.  For  this  reason,  we  restrict  this  discussion  to  GTO 
periods,  although  the  same  general  principles  apply  to  the  LASSII  model.  The  present  discussion  is  also 
restricted  to  attitude  characteristics  that  are  essential  in  the  modeling.  Other  features  will  be  discussed 
in  later  sections. 


3.1  A  TYPICAL  ORBIT 


A  plot  of  the  attitude  state  vectors  for  Orbit  31j,  typical  of  most  early  orbits  with  eclipse.  sliown  in 
Figure  5.  The  plot  is  from  the  attitude  quality  control  and  fitting  software  developed  by  Radex,  ...:d 
consists  of  four  panels.  The  top  panel  shows  the  attitude  source  fiag,  which  shows  tne  data  that  went 
into  the  calculation  of  the  attitude  at  each  point.  Tlie  values  taken  by  the  source  flag  are  the  same  as 
those  given  in  Table  1.  'Fhe  next  panel  gives  the  spin  axis  right  ascension,  the  next  panel  the  declination 
and  the  bottor.  pt.nel  the  spin  rate  in  rpm.  Concentrating  first  on  the  spin  axis  pointing  direction,  we 
can  see  that  tac  data  indicates  a  total  of  five  discontinuities  in  the  pointing  direction  as  well  as  a 
persistent  drift  \vhjch  seems  to  decrease  the  declination  by  about  5®  at  apogee. 

In  fact,  these  features  of  the  data  are  not  real.  Tliey  result  from  the  use  of  the  magnetic  field  model  and 
the  engineering  magnetometer  for  the  calcLLdio.-r  cf  .iltltude  in  periods  when  the  resolution  of  the 
magnetometer  and  the  acciirucy  of  the  assumed  model  held  are  not  sufficient  for  accurate  determination. 
One  source  of  the  discontinuities  is  jumping  of  the  attitude  solution  between  the  two  possible  solutions 
for  each  data  point.  Tliis  happens  when  the  first  guess  attitude,  which  is  substantially  more  accurate  than 
either  of  the  hvo  individu.ii  attitudes  ailculated  from  apogee  magnetic  field  hits,  falls  between  the  two. 
A  second  source  is  the  downlink  power  amplifier  which  causes  a  rather  large  current  loo|r  when  activated, 
and  induces  a  strong  spacecraft  magnetic  field.  In  the  absence  cf  complicating  factors,  we  would  expect 
the  spin  axis  pointing  direction  to  be  constant  between  attitude  adjustments.  There  are  in  fact  some 
compliciting  factors  for  CRRES,  to  be  discussed  later,  but  these  turn  out  to  be  too  small  to  be  included 
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Figure  5.  An  attitude  survey  plot  for  Orbit  313,  showing  the  calculated  right  ascension, 
declination  and  spin  rate  from  each  observation  set. 


in  the  model  on  a  routine  basis.  'Hie  most  accurate  attitude  in  Figure  5  is  that  given  by  tlie  data 
including  earth  horizon  hits  and  sun  hits  (source  code  7). 

A  second  feature  seen  in  Figure  3  is  a  rather  dramatic  increase  in  spin  rate  during  eclipse  (source  code 
4)  and  a  more  rapid  decrease  after  the  satellite  emerges.  Tliis  spin  rate  change  is  real  and  is  the  result, 
for  the  most  part,  of  cooling  and  heating  of  the  wire  antennae.  A  third  feature  which  is  not  so  readily 
apparent  in  this  orbit  is  a  change  in  spin  rate  or  spin  axis  pointing  direction  very  close  to  both  first  and 
second  perigee.  Because  of  the  interrelation  between  spin  axis  pointing  direction  and  spin  rate 
calculations,  it  is  difficult  to  differentiate  between  the  two  for  a  change  as  small  as  this.  ITie  change  is 
significant  enough,  though,  to  warrant  segmentation  of  the  near  earth  portions  of  the  orbit. 
Comparisons  with  measured  magnetic  field  have  shown  that  this  is  most  satisfactorily  treated  as  a  change 
in  spin  rate  rather  than  in  spin  axis  pointing  direction. 


3.2  CADP  FIT  RESUL'rS 


Figure  6  shows  an  attitude  quality  control  plot  for  the  attitude  model  returned  by  the  CADP  for  Orbit 
313.  This  plot  is  also  from  the  attitude  modeling  system  and  a  plot  such  as  this  is  produced  for  every 
orbit  processed.  The  top  three  panels  show  the  state  vector  values  of  a,  5  and  o  plotted  as  solid  lines 
and  the  model  values  as  dots.  The  bottom  panel  shows  the  phase  of  the  sun,  measured  in  the  spacecraft 
frame,  at  the  time  of  sun  hits.  Because  the  sun  sensor  is  on  the  -x  side  of  the  spacecraft,  this  value 
should  be  180°.  The  phase  is  calculated  as  follows.  First,  at  the  time  of  each  sun  hit,  the  instantaneous 
values  of  a,  5  and  u  are  calculated.  'Ihen,  the  spin  phase  is  found  from 

4*(iss)  ~  ^(^ss '  ^0^ 

The  90°  is  added  to  make  the  definition  of  the  CRRES  attitude  correspond  to  the  standard  rotation 
matrbc  below.  The  matrix  M  describing  the  instantaneous  pointing  directions  of  the  spacecraft  axes  in 
ECI  coordinates  is  given  by 


M  = 


sin6sinacos(i)-sinasin<|) 
-sin6cosasin(j)-sinacos(j) 
cos 6 cos  a 


sin6sinacos(t)+cosasin4) 

-sin5sinasin({)+cosacos{t) 

cos6sina 


-cos6cos({) 
cos6sin(j) 
sin  6 


(6) 


The  ECI  position  of  the  sun  s  is  interpolated  from  the  CRRES  ephemeris  then  the  sun  vector  in 
spaceaaft  coordinates  is  given  by 

Sj/c  =  M  s 

The  phase  is  then  calculated  from  (j)^  =  tan'’(Sy/Sjj). 
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Figure  6.  Quality  control  plot  of  the  altitude  model  compared  to  data  used  for  modeling 
for  the  original  CADP  fitting  routine. 


18 


As  can  be  seen,  the  CADP  attitude  model  has  several  undesirable  characteristics.  First,  the  values  of  a 
and  5  have  been  modeled  as  linear,  quadratic  and  even  quartic  polynomials  when  in  fact  they  should  be 
nearly  constant.  The  errors  in  the  calculated  attitude  away  from  perigee  have  even  been  magnified  by 
the  fitting.  One  would  prefer  to  use  only  the  reliable  data  for  modeling  of  a  and  6.  lliis  is  entirely 
possible  since  the  values  are  constant  throughout  an  orbit.  However,  this  doesn’t  fit  into  the  CADP 
modeling  scheme  since  an  independent  model  for  a  and  6  is  produced  for  each  segment  from  data  in 
that  segment  alone. 

Also,  during  the  eclipse  (labeled  BE)  and  the  post-eclipse  segment  (EE),  the  spin  rate  has  been  modeled 
as  a  constant.  We  can  see  from  the  phase  data  that  this  results  in  errors  of  as  much  as  one-hundred 
degrees  in  the  post-eclipse  segment.  (Phase  is  not  ailculated  for  this  plot  during  eclipse.)  Qeiirly,  a 
constant  spin  rate  will  not  do.  Finally,  looking  closely,  we  see  some  rather  large  discontinuities  in  phase 
at  segment  boundaries,  even  for  sunlit  inter\als.  'Fhis  am  be  seen  more  clearly  in  Figure  7,  which  show's 
the  angle  between  the  magnetic  field  measured  b>  the  science  magnetometer  and  the  model  field.  'Ilie 
measured  field  requires  conversion  to  ECI  in  order  to  calculate  this  angle  so  the  angular  deviation  am 
be  taken  to  be  the  combination  of  errors  in  ailibration,  magnetometer  alignment  and  iliree-axis  attitude. 
It  also  includes,  of  course,  the  difference  between  the  actual  magnetic  field  and  the  model.  We  will  ha\e 
more  to  say  about  this  comparison  later,  but  for  now,  it  is  apparent  that  discontinuities  of  10-20'^  are 
present  and  that  the  attitude  even  during  sunlit  segments  is  clearly  often  in  serious  error.  Discontinuities 
between  segments  of  more  than  one-quarter  degree  or  so  can  be  very  annoying,  especially  for  magnetic 
field  despinning  where  they  cause  low  frequency'  contamination  of  wave  spectrograms  and  jumps  in  the 
despun  magnetic  field. 


3.3  CONSTRAINTS 


Before  launching  into  the  solutions  adopted,  it  seems  a  good  idea  to  list  the  factors  limiting  the  modeling 
effort.  These  are  the  assumptions  that  have  been  made  in  attacking  the  improvement  of  the  model. 

Tht  aiiiiudt  must  be  modeled  based  on  pure  rotation.  Tliis  requirement  arises  for  the  most 
part  from  the  complexity  of  the  dynamical  motion  of  the  satellite.  It  was  never  considered  practical  to 
attempt  to  model  nutation  following  attitude  adjustments,  if  it  is  indeed  possible.  Likewise,  even  though 
it  was  known  for  some  time  that  the  wire  booms  would  interact  strongly  after  most  attitude  maneiners. 
it  was  alw'ays  considered  too  large  a  task  to  attempt  to  produce  anything  more  than  an  *.i\eraged’ 
attitude.  Fortunately,  modeling  the  attitude  in  this  way  allows  us  to  obtain  the  required  accuracy  most 
of  the  time. 

Orbits  must  be  processed  independently.  Because  of  variability  in  the  time  it  takes  for  tapes 
to  arrive  from  different  stations,  it  is  often  the  case  that  orbits  are  processed  out  of  sequence.  For  this 
reason,  w^e  can  never  count  on  being  able  to  use  data  from  a  previous  or  later  orbit  in  the  modeling  of 
attitude  for  any  particular  orbit. 

Attitude  must  be  processed  concurrently  with  Agency  Tape  Generation,  'iltis  limitation  w^as 
strongest  in  the  early  orbits,  when  methods  to  deal  with  non-nominal  behavior  of  the  attitude  were  still 
being  developed.  Even  later,  though,  when  for  e.xample  the  length  of  eclipse  increased  so  tliat  a 
quadratic  spin  rate  increase  provided  a  more  accurate  representation  than  a  linear  one,  some  orbits  were 
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Figure  7.  The  angle  between  a  model  magnetic  field  and  the  field  as  measured  by  the 
fluxgate  magnetometer.  Attitude  is  modeled  by  the  original  CADP  routine. 
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processed  while  algorithm  modifications  were  being  made  which  could  have  been  better  treated  using  the 
present  model.  This  impacts  the  user  in  that  he  should  bear  in  mind  that  all  attitude  is  not  equal  and 
that  in  some  instances  should  the  Agency  Tape  version  of  the  model  prove  inadequate,  a  better  model 
may  be  available. 

The  attitude  model  must  be  consistent  with  AGMOD.  In  dealing  with  the  spin  rate 
variations,  it  might  have  been  far  easier  to  redefine  the  spin  phase  used  for  modeling  the  attitude. 
However,  that  change  would  have  required  a  revision  of  the  software  subroutine  AGMOD  used  by  the 
agencies  to  evaluate  the  model.  While  this  would  not  in  all  likelihood  have  been  a  disaster,  and  would 
certainly  have  been  done  if  it  would  have  enhanced  the  accuracy  of  the  result,  it  seemed  simpler  to  leave 
the  structure  of  the  model  the  way  it  was.  Tlius,  AGMOD  has  only  a  single  version. 

The  attitude  model  must  be  based  on  data.  We  have  opted  to  avoid  heuristic  models  in  the 
CRRES  attitude,  at  least  so  far.  lliis  pertains  mostly  to  the  modeling  of  the  spin  rate  change  during 
eclipse,  since  data  there  is  limited.  One  reason  for  this  limitation  is  that  it  avoids  potential  problems 
with  operational  changes.  For  example,  during  several  orbits,  some  of  the  heaters  were  turned  on  during 
the  eclipse.  This  led  to  a  far  different  characteristic  in  the  spin  rate  change.  While  the  data  based  model 
used  may  not  have  been  as  accurate  here  as  in  unheated  orbits,  the  use  of  a  heuristic  model  would 
probably  been  wildly  incorrect  and  might  well  have  gone  unnoticed.  Another  more  practical  reason  was 
that  in  the  past  eclipse  season,  the  duration  of  eclipses  was  increasing  instead  of  decreasing  and,  short 
of  a  fully  analytical  model  of  spin  rate  change  which  is  likely  beyond  our  cipacity,  a  model  based  on 
previous  orbits  would  always  have  to  be  extrapolated. 

The  attitude  processing  must  not  slow  down  agency  tape  generation.  'Hiis  is  a  restriction 
in  that  it  limits  the  complexity  of  the  model.  Anything  we  do  must  be  successful  in  a  majority  of  the 
cases  with  a  minimum  of  intervention.  Since  attitude  determination  and  modeling  is  nearly  the  final  step 
in  tape  generation,  the  processing  is  done  under  some  pressure  of  time.  The  algorithms  developed  must 
perform  satisfactorily  on  the  vast  majority  of  orbits  so  as  not  to  cause  data  processing  to  fall  significantly 
behind. 
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4.  ATTITUDE  MODELING 


In  this  section,  the  algorithms  and  features  of  the  CRRES  attitude  model  will  be  described.  The 
modeling  itself  is  done  with  a  fully  interactive  software  package,  after  batch  processing  of  the  attitude 
state  vectors.  The  basic  features  of  this  package  have  been  discussed  before  [McNeil,  1990].  The  focus 
here  will  be  to  describe  the  methodology  and  specific  algorithms  that  were  developed  to  deal  with  the 
problems  outlined  in  Section  3. 


4.1  SPIN  AXIS  VECTOR 


As  previously  pointed  out,  the  only  reliable  spin  axis  data  comes  from  periods  of  the  orbit  with  both  sun 
and  earth  limb  contacts.  We  also  claimed  that  the  spin  axis  pointing  direction  was  constant  during  an 
orbit  in  the  absence  of  an  attitude  adjustment.  Tliis  is  not  precisely  true.  From  a  survey  of  a  large 
number  of  orbits  around  Orbit  250,  we  have  found  that  the  spin  axis  attitude  does  indeed  change  by 
something  on  the  order  of  0.05°  per  orbit.  There  are  two  reasons  for  ignoring  this  change.  First,  the 
change  is  too  small  to  be  detected  reliably  within  a  single  orbit  since  the  maximum  resolution  is  around 
0.5°  mainly  on  account  of  the  sun  sensor.  Since  we  are  required  to  process  each  orbit  independently, 
incorporation  of  this  small  change  into  the  modeling  would  be  impossible.  Also,  it  seems  likely  that  this 
change  is  due  to  a  gravity  gradient  and/or  atmospheric  drag  near  perigee.  If  so,  then  the  change  occurs 
over  a  rather  short  period  of  time  very  near  perigee  and  the  attitude  would  be  expected  to  be  quite 
constant  throughout  the  greater  part  of  the  orbit.  We  find  that  a  constant  spin  axis  pointing  direction 
is  adequate. 

Normal  Orbits.  A  typical  GTO  orbit  will  have  two  to  three  periods  of  earth  horizon  data.  Tlie 
number  and  duration  of  these  depends  on  the  orientation  of  perigee  with  respect  to  the  sun.  To  calculate 
the  averaged  values  of  a  and  ft  in  an  orbit,  a  global  fit  is  done  on  the  values  for  all  state  vectors  weighing 
very  heavily  the  data  from  sun  plus  earth  observation  sets.  These  values  are  used  in  all  segments  of  the 
orbit.  We  have  not,  to  date,  encountered  an  orbit  without  earth  data,  even  when  periods  of  data  are 
missing  due  to  missing  tapes.  Should  CRRES  go  to  partial  data  coverage  because  of,  for  example,  failure 
of  one  of  the  tape  recorders,  this  possibility  might  well  have  to  be  addressed.  One  solution  would  be  to 
relax  the  requirement  against  using  data  from  successive  orbits.  Since  the  change  in  pointing  direction 
is  quite  small  from  orbit  to  orbit,  values  from  a  previous  orbit  could  be  used  without  degrading  the 
accuracy  of  the  solution. 

Attitude  Adjustments.  In  orbits  with  attitude  adjusts. global  average  of  a  and  6  is  calculated 
with  data  both  before  and  after  the  adjustment.  For  these,  the  modeled  pointing  direction  changes 
discontinuously  at  the  beginning  of  the  first  adjustment  segment.  For  a  few  orbits  missing  first  perigee 
earth  hits  and  containing  attitude  adjustments  early  on  in  the  orbit,  the  lack  of  pre-adjust  data  presented 
a  problem.  For  these  cases,  sun  plus  magnetic  field  data  very  near  first  perigee  was  used  to  calculate  the 
pre-adjust  attitude.  'The  attitude  result  from  the  sun  and  magnetometer  within  a  couple  of  thousand 
kilometers  of  perigee  is  virtually  indistinguishable  from  that  ailculated  with  earth  data  anyway. 
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4.2  SPIN  RATC 


In  that  the  instantaneous  spin  phase  in  the  model  is  calculated  from  an  initial  phase  and  the  spin  rate, 
the  modeling  of  spin  rate  and  spin  phase  are  interchangeable.  As  noted  previously,  one  severe  limitation 
of  the  original  attitude  model  was  that  the  spin  rate  had  to  be  constant  within  a  segment.  Another 
problem,  not  readily  apparent  in  Figure  6,  was  that  the  spin  phase  zero  time  for  each  observation  set  was 
calculated  using  the  spin  axis  pointing  direction  for  that  observation  set.  We  have  seen  that  the  vast 
majority  of  these  are  in  error  by  several  degrees.  'Fhis  error  would  be  propagated  into  and  is  in  fact 
magnified  in  the  phase  zero  times  used  to  model  the  phase  if  they  were  not  corrected. 

Tlie  first  step  in  modeling  the  phase  is  to  use  the  global  averaged  right  ascension  and  declination 
determined  from  earth  sensor  periods  to  ailculate  a  new  set  of  phase  zero  times  from  the  time  of  each 
obser\'ation  set  containing  sun  data,  'flie  time  assigned  to  each  observation  set  is  also  the  time  of  the 
sun  hit.  Observation  sets  without  sun,  t.e.  during  eclipse,  are  treated  below,  llie  calculation  of  a  new 
set  of  phase  zero  times,  tQ  j  goes  precisely  as  described  in  Section  2.7  with  the  global  average  z  replacing 
the  observation  set  z  used  in  the  CADP. 

Next,  these  phase  zero  times  are  used  to  calculate  a  phase  variable  as  follows.  An  averaged  constant 
spin  rate  is  calculated  over  the  segment,  llien  the  phase  zero  times  are  used  to  calculate 

<l>i  =  "a(‘o,i  -  h) 

These  are  modulated  by  2n  and  corrected  should  the  value  pass  2jc  on  successive  points.  Since  the  (j)i’s 
represent  the  phase  plus  the  initial  phase  at  the  time  of  phase  zero,  if  a  plot  of  <}>j  against  tQ  j  is  constant, 
then  is  the  true  spin  rate  and 


4)q  =  -(j). 

Since  (|)  is  the  integral  of  ca,  the  determination  of  phase  from  phase  data  instead  of  from  spin  rate  data 
seems  likely  to  be  more  accurate  than  the  calculation  of  phase  from  a  set  of  successive  measured  spin 
rates.  From  the  variation  in  the  <j),  values,  constant  through  quadratic  correction  terms  to  the  spin  rate 
are  determined.  Since  o)  is  a  derivative,  this  means  that  the  model  allows  for  up  to  cubic  variation  in  the 
phase  during  a  segment. 

For  no  particular  re;».son  save  convenience,  the  determination  of  the  correction  terms  is  done  by  fitting 
4>,  with  the  same  Qiebyshev  polynomials  as  were  used  for  a  and  6.  We  then  am  obtain  the  phase  at 
any  point  X  from 


(t)(X)  =  6(0 +  Kq  +  7C  jX  +  “  1]  -  3X]  (7) 


where  the  re,  are  the  expansion  coefficients.  Evaluation  of  the  phase  by  the  attitude  model,  however, 
requires  the  expansion  of  the  spin  rate  as  follows 
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4)(t)  =(t)o  +6'c(a)o  +  WiA’  +  (i)2[2X^-l]) 


(8) 


where  T=t-tj.  When  evaluated  in  terms  of  x,  this  gives 

^2  ^3 

TT  T 

<|)(t)  =  4)o  +  6((i)q  -  tOj  +  (1)2)'^  +  6(20)^  -80)2) —  +  6-80)2 —  (9) 

A 

where  A=t2-tj.  As  may  be  seen,  it  would  have  been  substantially  easier  to  expand  the  phase,  modulated 
by.the  average  spin  rate,  in  Qiebyshev  coefficients  in  the  modeling.  Since  the  agency  module  expands 
the  spin  rate,  however,  we  must  determine  the  coefficients  o,  from  the  measured  coefficients  rCj.  Before 
doing  so,  though,  we  must  note  one  complication.  Because  the  quantity  fit  by  the  Chebyshev  expansion 
is  actually  the  value  required  to  obtain  zero  phase  at  each  fit  time,  the  true  phase  is  given  by 

(l)(;0  =  6o)^t  -  TCo  -  It  j X  -  •n:2[2A’^  - 1]  -  -  3X]  (10) 


First,  evaluating  the  Eq(lO)  at  t=0  and  X=-l,  we  find 

^0  “  “^0  “  ^2  ^3 


(11) 


Next,  we  evaluate  Eq(7)  and  Eq(9)  at  t=0.25A  and  X=-0.5  to  give 


0)0  0)i  0)2^ 


3o)^A  -Ttj  +3Tt2-4it3  =  6A( —  - —  - — ) 

2  4  4 


(12) 


In  the  same  way,  evaluation  at  t=A  and  X=1  gives 

6o)^A  -27rj  -27:3  =  6A(o)o+o)j+o)2) 


(13) 


and  evaluation  at  t=A/2  and  X=0  gives 


0)o  0)2^ 


3o)^A  -TCj  +2Tt2-TC3  =  6A( —  - — ) 


Taking  4  times  Eq(12)  gives 

12o)^A  -4it,  +121:2-16713  =  6A(2o)o-o)j  -0)3) 


(14) 


(15) 
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Adding  this  to  Eq(13)  we  find 


(Oq  = 


18(i)^A  -GtCj  -4it2-l8it3 
_ 


(16) 


Next,  from  Eq(14)  we  find 


“2  =  “0- 


6(0^ A  -271j  +4it2-27t3 

6A 


(17) 


and  from  Eq(13)  we  find 


(i)j  =  -(j)q-(»)2  + 


6(»)^A  -2Ttj  -2it3 

6A 


(18) 


4.3  ECLIPSE  MODEL 


During  the  eclipse,  the  only  information  of  any  practical  use  is  the  spin  rate  measured  by  the  time 
between  the  zero  crossings  of  the  x-axis  and  the  magnetic  field.  To  model  the  spin  rate  increase  during 
eclipse,  we  perform  a  fit  to  ihe  spin  rate  using  the  Chebyshev  expansion.  We  then  use  characteristics 
of  this  fit  to  model  the  phase. 

Assume  that  a  quadratic  fit  of  w  during  the  eclipse  results  in  coefficients  Oq,  and  Oj.  We  use  the  first 
and  second  derivatives  of  the  measured  o)(t)  for  the  first  and  second  derivatives  of  the  spin  rate  in  the 
model.  These,  as  with  the  previous  correction  terms,  require  modification  because  of  the  way  in  which 
phase  is  evaluated  by  AGMOD.  llie  resulting  coefficients  are 

<i)j  =  o,  -  4a^6 

©2  =  02^3 

We  match  the  phase  at  the  start  of  the  eclipse  segment  to  the  phase  at  the  end  of  the  previous  segment. 
We  modify  the  constant  term  in  the  spin  rate  according  to 

“o  “  “c  ■  *^0  ■ 

where  is  a  correction  term  to  make  the  phase  at  the  end  of  the  eclipse  segment  match  the  phase  at 
the  beginning  of  the  next  segment,  calculated  as  follows.  If  a  is  the  number  of  ‘excess’  revolutions  in 
the  segment  calculated  from  the  difference  between 
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and  the  phase  at  the  beginning  of  the  next  segment,  then 


<i)j.  =  60a/(t2  -  tj) 

4 

with  CD,  as  usual,  in  rpm. 

Finally,  because  the  segment  following  the  eclipse  is  very  short  and  because  wire  boom  motion,  to  be 
discussed  later,  often  confounds  the  calculations  somewhat  in  this  area,  the  phase  at  the  end  of  the  post¬ 
eclipse  segment  (EE)  is  matched  to  that  at  the  beginning  of  the  next  segment  by  modification  of  the 
constant  term  in  the  spin  rate.  ITiis  is  done  in  the  same  way  as  in  the  previous  calculation. 

We  should  mention  one  further  point  concerning  the  transformation  of  spin  rate  to  phase.  Because  the 
spin  rate  expansions  in  Sections  4.2  and  4.3  consist  of  terms  which  wiU  give  the  correct  of  the  phase 
instead  of  the  spin  rate,  evaluation  of  the  Chebyshev  expansion  with  these  coefficients  will  not  give  the 
correct  spin  rate.  The  expression  for  the  instantaneous  spin  rate  using  the  spin  rate  coefficients  is 

(i)(A)  =  (i)q  -  <i)j  -1-  (Dj  +  (4(i)j  -  16ci)2)A  -b  24ci)2A^ 

where 

A  =  (t-ti)/(t2-ti) 

'I'his  is  of  no  practical  use  unless  one  wants  to  calculate  the  instantaneous  spin  rate  given  by  the  attitude 
model,  as  we  do  in  the  survey  plots. 


4.4  SEGMENTATION 


Because  of  the  incorporation  of  the  spin  rate  correction  terms  based  on  phase  fitting  into  the  modeled 
spin  rate,  it  is  possible  to  increase  the  duration  of  segments,  at  least  in  those  periods  for  which  a 
quadratic  spin  rate  is  a  good  approximation.  'Fhis  has  the  advantage  of  reducing  the  number  of 
discontinuities  in  an  orbit.  The  modeling  software  automatically  segments  normal  orbits  with  a  single 
segment  up  to  eclipse  time  (if  there  is  one),  and  two  to  three  segments  after  eclipse.  This  is  to  provide 
an  adequate  representation  of  the  spin  rate  decrease  after  emergence  from  eclii)se.  In  every  orbit  small 
segments  are  placed  before  and  after  perigee,  unless  this  time  comes  during  eclipse.  This  is  to  allow  for 
more  accurate  modeling  of  what  appears  to  be  a  spin  rate  change  at  perigee.  Also,  segments  are  placed 
interactively  where  mission  events  or  anomalies  occur.  A  survey  of  these  anomalies  is  given  in  Section 
5. 

We  should  emphasize  that  the  segmentation  varies  depending  on  the  characteristics  of  a  particular  orbit. 
Several  orbits  prior  to  Orbit  300  were  processed  with  two  post-eclipse  segments  and  a  linear  variation 
of  the  spin  rate  in  eclipse.  This  does  not  mean  that  they  are  substantially  less  accurate,  hov>cver,  since 
the  shorter  eclipses  gave  smaller  spin  rate  changes  and  required  fewer  segments. 


26 


By  orbit  425,  eclipse  had  disappeared  and  orbits  without  anomalies  of  other  sorts  were  treated  with  two 
one-hour  segments  at  either  end  and  a  single  segment  for  the  remainder.  When  eclipse  aime  back  again, 
around  orbit  595,  the  old  eclipse  model  was  used  for  a  while.  Then,  post-eclipse  segments  were  modified 
or  removed  when  eclipse  slipped  nearer  second  perigee.  By  orbit  800,  subtle  changes  in  the  spin  rate 
behavior  forced  the  addition  of  a  quadratic  term  to  the  spin  rate  in  sunlit  periods. 

'fhis  rather  rough  history  of  the  chosen  segmentation  is  intended  to  show  that  the  attitude  model  is 
anything  but  static,  the  characteristics  dependant  to  a  large  extent  on  the  judgment  of  the  attitude  analyst, 
within  the  constraints  of  the  specified  tolerable  error  level  and  other  constraints  enumerated  previously. 
Schemes  are  modified  when  it  becomes  apparent  that  the  current  method  is  becoming  inadequate  and 
when  reasonable  changes  will  lead  to  a  significantly  better  result. 


4.5  RESULTS 


Figure  8  shows  an  attitude  quality  control  plot  for  the  model  calculated  as  described  above.  The 
corresponding  plot  of  the  angular  deviation  of  the  measured  and  model  field  is  shown  in  Figure  9.  As 
should  be  apparent  from  Figure  8,  the  attitude  model  reproduces  the  data  quite  well,  with  the  exception 
of  a  and  5  at  off-perigee  times.  As  we  know,  the  calculated  values  of  a  and  5  are  incorrect  except  when 
there  are  earth  hits  and  very  near  perigee,  where  the  modeled  values  intercept  the  data  as  they  should. 
The  modeled  spin  rate  follows  the  calculated  one  quite  well  and  the  phase  calculation  indicates  that  the 
model  phase  is  consistent  with  the  data  to  within  one-half  degree  at  most.  The  19  minute  oscillation  in 
the  values  of  the  sun  phase,  given  in  the  bottom  panel,  arise  from  the  oscillations  of  the  wire  booms. 
Since  the  sun  is  high  in  the  spacecraft  frame,  the  actual  motion  of  the  spin  axis  is  substantially  less  than 
would  be  indicated  by  the  magnitude  of  this  oscillation. 

Turning  to  the  magnetic  field  plot,  shown  in  Figure  9,  we  see  that  the  angular  deviation  between 
measured  and  model  field  is  less  than  about  2°  throughout  the  orbit.  Tliis  orbit  is  fairly  typical  as  far 
as  attitude  determination  is  concerned  but  is  somewhat  unique  in  that  magnetic  field  is  unusually  quiet, 
leading  to  good  agreement  between  the  measurements  and  the  model  field,  which  is  IGRF85  extrapolated 
to  date  plus  01son-Pfitzer77  Quiet  External  Field  Model.  As  mentioned,  the  angular  deviation  between 
measured  and  model  field  has  components  arising  from  magnetometer  calibration  and  alignment,  true 
discrepancy  between  the  model  and  actual  field  and  attitude  errors.  Had  we  chosen  a  more  disturbed 
period,  we  might  see  something  like  Figure  10.  'lliere,  the  agreement  is  still  less  than  2°  near  perigee, 
but  deviates  by  several  degrees  once  the  external  magnetic  field  comes  into  play,  lliese  deviations 
underscore  the  fact  that  magnetic  field  based  attitude  calculation  is  of  little  use  more  than  a  few  thousand 
kilometers  from  the  earth.  The  attitude  modeling  does  not  suffer  from  this  limitation  because  we  have 
chosen  to  rely  on  sun  and  earth  sensor  data  for  attitude  ailculafion. 

In  spite  of  the  large  deviations  caused  by  the  inadequacy  of  the  model  field,  we  can  still  use  this  plot  to 
evaluate  the  soundness  of  the  attitude  throughout  the  entire  orbit.  First,  we  note  that  the  spin  axis 
pointing  direction  is  the  same  throughout  the  orbit  and  the  phase,  based  on  timing  of  sun  hits,  should 
be  equally  accurate  in  all  sunlit  segments.  Tlius,  the  attitude  points  near  apogee  should  be  of  the  same 
level  of  precision  as  those  at  perigee.  Next,  we  note  that  the  attitude  is  continuous  throughout  a  segment 
but  generally  discontinuous  at  the  boundaries.  Because  of  the  general  trend  of  the  deviation  from 
segment  2  until  10:20  when  the  waves  are  encountered,  we  can  be  safe  in  mentally  subtracting  out  what 
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CRRES  POP  ATTITUDE  FROM  CflF03l30 
ORBIT  313  DAY  335  YEAR  1990  UTO  33780.00 


RASC 


1991-08-05 


TIME  FROM  START  OF  ORBIT  (SEC) 


Figure  8.  Same  as  Figure  6  but  using  the  model  generated  by  the  present  attitude 
modeling  software. 
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)99l-0e-05  16^322  V.2 

CRRES  Attitude  CXjdity  Control  Plot  313  335/90  Dec  1 


R(gr.)  2.18  4.05  5.24  5.94  6.25  6.19  5  75  4.90  3.51  1.39 

Lot(mag)  7.3  17.3  17.9  16.6  14.6  12.2  9.3  6.0  1.1  -8.6 

ULT  21.59  0:24  1:28  2.11  2:47  3.22  3.59  4:46  6:02  10:30 

2.15  4.47  5.95  6.70  6.90  6.65  6.03  5.02  3.51  1.42 


Figure  9.  Same  as  Figure  5,  but  witii  the  attitude  model  generated  by  the  current  attitude 
modeling  software. 


29 


Figure  10.  Angular  deviation  between  ineabured  and  model  field  during  a  disturbed 


period. 
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is  likely  a  natural  deviation.  Another  factor  reinforcing  this  is  that  the  maxiimiin  errors  in  a  segment 
should  come  at  the  segment  boundaries.  Since  the  segments  are  fit  independently,  the  level  of 
discontinuity  at  the  boundaries  is  a  measure  of  the  maximum  error  expected.  Here  discontinuities  are 
all  less  than  one  degree.  The  good  match  between  the  segment  beginning  at  13:00  and  the  one  at  15:30 
indicates  as  well  that  the  attitude  during  these  periods  is  probably  no  worse  than  it  is  at  perigee.  During 
the  BE  and  EE  segments,  we  can  once  again  rely  on  the  actual  value  of  the  deviation  since  model  and 
measured  field  are  expected  to  agree  well  here. 


5.  ERROR  ESTIMATES 


In  this  section,  estimates  are  given  for  the  expected  level  of  accuracy  of  the  CRRES  attitude  model. 
Accuracy  is  judged  mainly  from  plots  like  Figures  9  and  10  of  the  angle  between  the  magnetic  field 
measured  by  the  science  magnetometer  and  the  IGRF85  model  field.  Tlie  calibration  and  alignment  of 
the  magnetometer  is  therefore  crucial  in  this  comparison.  The  elevations  and  relative  phase  angles  of 
the  magnetometer  axes  are  measured  routinely  during  on-orbit  calibration.  However,  the  absolute  phase 
of  the  instrument  is  somewhat  in  question  due  to  mis-deployment  of  the  magnetometer  boom,  'llie  value 
chosen  here,  177.4°,  was  calculated  by  minimization  of  the  angular  deviation  of  the  measured  and  model 
field  in  the  first  hour  after  perigee  over  the  course  of  many  orbits.  Variation  of  this  angle  would  more 
or  less  lead  to  the  addition  of  an  equivalent  angular  error  to  these  plots.  Tlie  important  point  to  bear 
in  mind  is  that  the  magnetometer  calibrations  and  alignments  are  the  same  for  each  plot.  Barring  some 
systematic  error,  t.g.  the  relative  timing  of  attitude  and  magnetometer  signals  (this  has  in  fact  been  ruled 
out  on  other  grounds)  these  plots  should  provide  reliable  error  estimates. 


5.1  NOMINAL  ORBrrS 


Tlie  vast  majority  of  the  magnetic  field  plots  show  deviations  like  those  of  Figure  9  or  10.  From  this, 
we  conclude  that  the  accuracy  of  the  vast  majority  of  CRRES  attitude  is  around  1-2°.  Exceptions  to  this 
arise  most  commonly  following  attitude  adjustments.  For  these  orbits  and  a  few  following  them,  errors 
in  certain  regions  exceed  this  level.  These  will  be  discussed  in  detail  below. 


5.2  ATHTODE  ADJUST  ORBI'IS 


Attitude  adjustments  are  made  to  keep  the  satellite  spin  axis  between  5  and  15°  of  sun  pointing. 
Problems  arise  following  the  adjustment  because  of  two  types  of  motion.  Nutation  of  the  spacecraft 
makes  it  difficult  in  some  cases  to  accurately  calculate  the  phase  after  the  maneuver.  Also,  a  very  strong 
interaction  of  the  spacecraft  with  the  wire  antennae  causes  a  long  period  oscillation  in  the  ailculated 
phase.  This  makes  for  difficult  modeling  of  periods  when  the  spacecraft  experiences  actual  spin  rate 
changes. 
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Since  the  attitude  model  for  CRRES  is  an  ‘averaged’  model,  meaning  that  it  does  not  attempt  to  include 
short  duration  periodic  motion  such  as  nutation,  we  must  first  examine  the  magnitude  of  such  motion 
following  attitude  adjustment  before  evaluation  of  the  ‘error’  in  the  model  itself.  Figure  11  shows  the 
measured  angle  between  the  sun  and  the  spacecraft  spin  axis  during  and  after  an  attitude  adjustment. 

We  see  nutation  of  as  much  as  10®  immediately  after  the  adjustment,  about  three-quarters  the  way 
through  Orbit  619.  By  the  beginning  of  Orbit  620,  this  has  damped  to  approximately  5®  and  is  about 
3®  by  the  beginning  of  Orbit  621.  In  what  follows,  then,  we  must  keep  in  mind  that  the  errors  shown 
in  the  magnetic  field  plots  represent  a  combination  of  errors  in  the  averaged  attitude,  which  is  the 
attitude  modeled,  and  true  short  term  oscillation  of  the  spacecraft  attitude. 

Figures  12  through  14  show  a  typical  attitude  adjustment  and  ^vo  orbits  following  it.  The  adjustment 
itself  takes  place  in  Orbit  802.  TTiis  orbit  is  further  complicated  by  three  clock  jumps  (TD),  one  about 
two  hours  after  the  maneuver.  From  the  magnitude  of  the  discontinuity  at  the  beginning  of  the 
maneuver  segment  (AA)  and  from  the  general  behavior  of  the  magnetic  field  before  the  maneuver  it 
would  seem  that  the  attitude  immediately  after  the  adjust  is  not  at  all  bad,  perhaps  in  error  by  5®  or  so. 
At  the  TD  segment,  though,  there  is  a  discontinuity  of  almost  20®.  Qock  jumps,  discussed  more  fiilly 
in  Section  V,  might  be  expected  to  give  rise  to  no  discontinuity  when  properly  processed  and  compared 
to  a  geophysical  requiring  attitude  for  the  comparison.  However,  the  wire  boom  motion  immediately 
following  the  maneuver  has  probably  led  to  improper  fitting  of  the  phase  near  the  segment  boundaries, 
lliis  happens  when  the  polynomial  fitting  algorithm  mistakes  a  partial  oscillation  near  segment  ends  for 
a  true  spin  rate  change.  The  smaller  the  segment  the  more  intense  this  effect  will  be  since  partial  cycles 
of  the  19  minute  oscillation  will  constitute  a  greater  portion  of  the  whole,  llie  same  effect  can  be  seen 
in  the  first  TD  segment,  which  has  discontinuities  in  the  2-3®  range.  Tliere,  the  segment  was  extremely 
short  due  to  a  second  clock  jump,  llie  oscillations  present  there  were  remnants  of  another  attitude 
adjust  some  two  orbits  prior  to  this  one.  Again  from  the  looks  of  the  overall  variation  of  the  field,  we 
must  admit  to  an  approximately  10®  error  near  the  beginning  of  the  third  TD  segment.  We  should 
emphasize  that  most  attitude  adjusts  do  not  contain  discontinuities  of  this  magnitude.  In  the  usual  case, 
a  segment  would  not  be  put  in  so  soon  after  the  maneuver. 

Moving  on  to  the  BE  segment,  we  can  see  again  the  ‘end  effects’  of  the  wire  boom  motion  leading  to 
errors  in  the  phase  near  the  segment  boundaries.  As  discussed  in  the  modeling  section,  the  phase  at  the 
beginning  of  the  BE  segment  is  matched  to  the  end  of  the  previous  segment  then  the  spin  rate  is 
adjusted  to  ir..atch  the  phase  at  the  beginning  of  the  EE  segment.  Tins  means  tliat  the  error  from  fitting 
the  third  TD  segment  is  carried  over  to  the  beginning  of  the  BE  segment.  Since  the  same  ‘end  effects’ 
are  present  in  the  EE  segment,  especially  severe  here  beaiuse  the  segment  is  so  short,  the  beginning  and 
end  phase  there  is  also  in  error  by  as  much  as  10®.  Of  this,  perhaps  one-half  is  actual  motion  of  the 
spacecraft  that  deviates  from  the  average  model  sought.  l*his  case  is  one  of  the  most  severe,  due  to  the 
very  short  EE  segment  and  due  to  the  presence  of  the  clock  jump.  However,  it  serves  the  purpose  in 
demonstrating  that  dynamical  motion  of  the  spacecraft  make  attitude  modeling  exlremel)  difficult  after 
an  attitude  adjustment.  We  cannot  increase  the  segment  durations  due  to  ecli[)ses,  clock  jumps  and 
simply  the  end  of  the  orbit.  Nor,  however,  can  we  accurately  ascertain  the  spin  rate  changes  because  of 
the  large  amplitude  oscillations  that  appear  due  to  wire  boom  interactions.  Fortunately,  this  interaction 
dampens  quite  rapidly,  as  we  will  see  shortly. 
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CRRES  Attitude  Qudity  Control  Rot  820  178/91  .Jun  27 


183  5.85  5.15  5S5  a38  6.<-;  5.15  550  4.39  2.67 


Lot(mog)  “”S  "2.1 

^gj  ’2-.20  1527  15:36  1720  1755  ^27  1359  1956  2023  22:09 

i/p,^  LSI  4.01  5.46  655  681  6.93  5.70  6.02  4  72  2.67 


Figure  12.  Angular  deviation  plot  for  an  orbit  containing  an  attitude  adjustment. 
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CRRES  Attitude  Quality  Control  Plot  821 179/91  Jun  28 
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Figure  13.  Angular  deviation  plot  for  one  orbit  after  an  attitude  adjustment. 
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CRRES  Attitude  Quolity  Control  Pbt  822  179/91  Jun  28 
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Figure  14.  Angular  deviation  plot  two  orbits  after  an  attitude  adjustment. 
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Figure  13  shows  the  orbit  following  the  attitude  adjustment  under  consideration.  First,  note  that  the 
magnitude  of  the  wire  boom  motion,  seen  most  clearly  around  11:30,  is  still  around  4®.  Next,  reaill  that 
in  the  first  segment,  we  are  trying  to  model  a  relatively  rapid  spin  rate  decrease  due  to  emergence  from 
eclipse  just  prior  to  perigee.  We  see  that  in  the  first  segment,  the  wire  boOm  motion  has  confounded 
the  fitting  of  the  phase  somewhat,  leading  to  an  initial  deviation  of  10®.  We  remember  again,  though, 
that  the  spacecraft  attitude  is  experiencing  a  complex  19  minute  oscillation  of  perhaps  5®  at  this  point. 
This  rapidly  drops  below  the  5®  level,  though,  at  least  on  average  when  the  actual  motion  of  the 
spacecraft  is  discounted.  In  segment  2,  we  see  errors  of  approximately  3®  and  those  in  segment  3  are 
probably  no  larger.  The  eclipse  as  well  shows  no  particular  changes  indicating  increased  error.  In  the 
final  EE  segment,  however,  we  see  again  that  wire  boom  interaction  has  been  translated  into  spin  rate 
changes  in  the  model  and  deviations  of  perhaps  as  much  as  5®  are  present  in  this  segment.  'Ihis  is  due 
mostly  to  the  ‘end  effects’  since  actual  nonmominal  motion  is  down  to  about  the  1®  level  at  this  point. 

Finally,  we  show  in  Figure  14  the  next  orbit,  beginning  about  one  and  one-half  orbits  after  the 
adjustment.  We  see  that  the  errors  in  the  first  segment  are  below  the  2®  level  and  those  in  the  remainder 
are  even  smaller,  near  1®.  Unfortunately,  the  EE  segment  in  this  plot  is  empty  due  to  a  switch  to 
LASSII  mode  and  consequent  loss  of  the  science  magnetometer  data.  Other  orbits,  however,  show  that 
this  error  is  less  than  2®  as  well. 


5.3  OTHER  ANOMALIES 


Errors  arising  in  the  calculation  and  modeling  of  the  attitude  during  other  periods  of  non-nominal 
behavior  are  not  particularly  severe.  Tliis  is  for  the  most  part  because  the  wire  boom  interaction 
resulting  firom  such  things  as  spin  rate  adjustments  or  canister  ejections  is  relatively  high  frequency  and 
does  not  interfere  with  spin  rate  modeling.  The  modeling  produces  an  averaged  attitude  as  desired 
without  contamination  from  boom  interactions.  The  only  aspect  that  might  be  of  concern  with  other 
types  of  anomalies  is  the  precise  timing  of  these  events.  At  the  time  of  attitude  processing,  there  are 
some  problems  inherent  in  the  determination  of  the  precise  times  of  attitude  events.  Pass  plans  are  often 
incomplete  and  mission  event  reports  are  often  not  available.  Consequently,  in  order  to  keep  pace  with 
processing,  the  times  for  all  unusual  attitude  events  and  the  segmentation  of  the  attitude  model 
accordingly  are  obtained  firom  inspection  of  the  attitude  data  itself.  Tliis  pertains  to  attitude  adjustments, 
spin  rate  changes  and  canister  releases,  as  well  as  the  timing  discontinuities  discussed  above.  ITie  most 
sensitive  measure  of  the  first  three  of  these  is  the  calculated  spin  rate,  available  only  at  the  spin 
frequenry,  since  it  is  extremely  sensitive  to  variations  in  phase  and  spin  axis  direction,  lliis  means  that 
the  segmentation  of  orbits  around  events  is  accurate  only  to  one  minute  or  so.  As  a  consequence, 
maneuvers  that  cause  a  change  in  the  spin  rate  may  result  in  attitude  models  with  phase  errors  of  a  few 
degrees  within  a  minute  of  the  exact  maneuver  time.  In  special  cases  where  the  data  immediately 
surrounding  an  attitude  maneuver  is  critical,  more  precise  segment  boundaries  can  be  incorporated  into 
the  model  and  a  revised  set  of  model  coefficients  issued  on  an  individual  bases,  lliis  has  been  done,  for 
example,  for  photometer  analysis  following  the  canister  releases. 
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5.4  SUMMARY 


The  above  examples  present  what  we  believe  to  be  an  honest  assessment  of  the  error  levels  to  be 
expected  in  the  CRRES  attitude  model.  As  a  general  rule  of  thumb,  one  can  expect  accuraqr  of  better 
than  2®  in  periods  more  than  two  orbits  after  an  attitude  adjustment.  An  attitude  adjust  orbit  may  have 
errors  of  10°  or  so  at  times  after  the  adjustment.  In  the  following  orbits,  segments  less  than  about  one 
hour  in  duration  may  also  exhibit  increased  error  levels  as  could  the  eclipse  segment  if  present.  Longer 
segments  should  for  the  most  part  be  accurate  to  2°  or  so.  Of  course,  there  are  exceptions  to  this  due 
to  special  events  or  orbital  configurations,  including  periods  of  data  dropout.  In  the  consideration  of 
accuracy,  it  is  also  important  to  keep  in  mind  that  we  are  discussing  the  accuracy  of  an  ‘average’  attitude. 
Angular  deviations  due  to  nutation,  in-plane  motion,  rocking  should  be  added  to  these  estimates. 
Approximate  descriptions  of  the  types  of  motion  to  be  expected  from  several  maneuvers  is  given  in  the 
next  section. 

In  support  of  this,  we  have  performed  a  survey  of  about  100  orbits,  calculating  the  RMS  average  and 
maximum  deviation  of  the  measured  field  from  the  IGRF  model  during  the  first  45  minutes  after  first 
perigee.  The  results  are  shown  in  Figure  15. 

All  of  the  orbits  around  625  and  690  with  RMS  deviations  exceeding  2°  are  immediately  following 
attitude  adjustments.  In  the  625  series,  the  adjustments  took  place  approximately  two-thirds  of  the  way 
into  the  orbit,  meaning  that  the  45  minutes  following  first  perigee  of  the  following  orbit  was  only  about 
one-third  orbit  after  the  adjustment.  As  we  have  seen,  deviation  of  the  spacecraft  attitude  from  pure 
rotation  can  account  for  about  5°  or  more  of  this  error,  so  even  a  perfect  model  would  still  fall  short  by 
this  amount. 

In  the  690  series,  the  attitude  adjustment  took  place  in  the  first  one-third  of  the  orbit  previous  to  the 
perigee  for  which  the  determination  in  Figure  15  was  made.  For  these,  the  lower  error  level  shows  the 
decay  of  the  wire  boom  motion  coupled  with  increased  ability  of  the  model  to  accurately  calculate  spin 
rate  changes. 

In  considering  Figure  15,  we  should  also  remember  that  the  determination  was  made,  by  necessity,  near 
first  perigee,  which  is  represented  by  a  short  segment.  Error  levels  in  the  much  longer  central  segment 
of  these  orbits  would  be  comparable  to  those  in  the  perigee  segments  of  orbits  for  which  wire  boom 
motion  was  not  a  significant  factor,  since  the  fitting  of  the  longer  segments  averages  over  the  motion  of 
the  wire  booms. 
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RnguLor  OevLaiLon 


Rnqular  DevLatLon  of  FLuxgale  Measured  FLeLd  and 
IGRF  during  the  First  -45  Munutes  of  Orbits  600~700 


Figure  15.  Deviation  of  the  measured  from  model  field  during  the  first  45  minutes  of 
Orbits  600-700. 


39 


6.  ANOMALIES 


This  section  describes  several  of  the  non-nominal  aspects  of  the  CRRES  attitude,  arising  either  from 
telemetry  problems  or  spacecraft  dynamics.  The  purpose  here  is  to  describe  generally  the  motion 
resulting  in  the  spacecraft.  Also  included  is  a  more  in  depth  discussion  of  clock  jumps. 


6.1  CLOCK  JUMPS 


From  time  to  time,  with  frequency  that  depends  do  some  extent  on  magnetic  activity,  the  spacecraft 
experiences  ‘clock  jumps’  in  which  the  Vehicle  Time  Correlation  Word  (VTCW)  loses  synchronization 
with  the  telemetry  formatting  system.  This  results  in  a  period  of  more  or  less  meaningless  telemetry  until 
the  clock  is  automatically  reset,  within  a  few  seconds  of  the  error.  This  has  two  important  effects  on 
attitude  processing.  First,  since  there  is  often  a  real  time  discrepancy  between  the  reset  value  of  the  clock 
and  the  passage  of  time,  there  appears  to  be  a  discontinuity  in  time.  For  example,  sun  sensor  data  timed 
from  the  VTCW  of  each  MF  might  indicate  the  passage  of  40  seconds  between  sun  hits  for  a  single  pair 
of  points  on  either  side  of  the  clock  jump.  All  other  pairs  would  be  close  to  the  30  seconds  expected 
at  2  rpm.  This  is  handled  in  the  attitude  model  by  placing  a  segment  near  the  clock  jump.  In  spite  of 
the  fact  that  the  spin  phase  will  be  discontinuous  at  the  segment  boundary,  this  method  allows  the 
attitude  model  to  match  the  rest  of  the  science  data  except  for  a  maximum  of  30  seconds  before  and  after 
the  jump.  The  uncertainty  region  of  maximum  60  seconds  arises  because  the  time  of  the  clock  jump  is 
determined,  as  for  all  mission  events  effecting  attitude,  from  the  attitude  data  itself.  The  time  resolution 
of  attitude  based  on  a  full  set  of  observations  is  one  spin  period.  Qock  jumps  are  seen  as  otherwise 
inexplicable  discontinuities  in  phase.  Although  more  precise  times  for  these  jumps  are  sometimes 
available  from  MFF  processing,  the  data  at  the  jump  is  invariably  bad  anyway  and  reliance  on  attitude 
data  alone  streamlines  the  processing  greatly. 

A  second  effect  of  these  jumps  is  that  the  absolute  UT  assigned  to  a  particular  MF  is  less  accurate  than 
that  of  the  time  correlations  from  which  it  is  derived.  Time  correlation  factors  match  a  particular  VCTW 
value  to  actual  UT.  Thus,  if  a  clock  jump  happens  between  two  measured  correlation  factors,  the  slope 
of  the  VCTW  to  UT  conversion  becomes  incorrect,  leading  to  incorrect  UT  values.  Again,  the  most 
sensible  way  to  deal  with  this  appears  to  be  to  make  the  times  of  the  attitude  model  match  the  timing 
of  the  rest  of  the  data,  allowing  that  the  actual  time  may  be  up  to  a  few  seconds  incorrect  in  some  cases. 
One  additional  complication  arises  from  the  need  to  relate  measured  data  to  ephemeris  data,  which  is 
timed  independently.  In  only  one  case  so  far  have  we  encountered  time  correlations  which  were  bad 
enough  to  cause  noticeable  discrepancy  beftveen  measurements  and  ephemeris  models.  'Fhese  few  orbits 
were  reprocessed  to  satisfaction  with  revised  time  correlation  factors. 
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6.2  ATTITUDE  ADJUSTMENTS 


Attitude  adjustments  are  carried  out  frequently  to  keep  the  spin  axis  within  5°  and  15°  of  sun  pointing. 
They  lead  to  strong  nutation  which  can  perhaps  15°  in  magnitude  initially.  The  dominant  frequencies 
are  approximately  once  and  twice  the  spin  period.  The  nutation  dampers  quench  the  nutation  down  to 
about  1°  three  hours  after  the  maneuver  and  to  about  0.1°  six  hours  later.  Attitude  adjusts  (and  all  other 
maneuvers)  are  modeled  by  an  ‘averaged’  attitude  in  the  CRRES  software.  No  attempt  is  made  to  model 
the  nutation  itself. 

'Ilie  energy  transferred  to  the  wire  booms  by  attitude  adjustments  auise  a  strong  oscillation  in  the  spin 
jixis  with  a  period  of  approximately  19  minutes.  'Hie  motion  consists  mainly  of  a  sort  of  rocking  of  the 
spin  axis.  TTie  magnitude  of  the  displacement  is  near  the  10°  level  up  to  one-half  orbit  or  so  following 
the  adjustment.  One  orbit  later  it  is  down  to  around  3°.  It  is  one-tenth  degree  three  orbits  or  so  later 
and  still  visible  in  magnetic  field  data  six  orbits  after  the  adjustment.  This  motion  aiuses  some  unique 
problems  in  attitude  modeling  beaiuse  of  its  long  period.  Spin  rate  changes  that  occur  on  timescales  of 
less  than  one-hour  or  so  are  quite  difficult  to  model  with  high  accurac7  when  the  wire  boom  motion  is 
intense.  This  is  because  the  spin  phase  ailculated  from  sun  sensor  data  will  show  a  low  frequency 
oscillation,  difficult  to  separate  from  the  actual  spin  rate  change.  Figure  11  shows  the  magnitude  of  this 
interaction  quite  well. 


6.3  SPIN  RATE  CHANGES 


Spin  rate  changes  are  made  whenever  the  spin  rate  drifts  too  far  from  the  nominal  2  rpm,  required  by 
several  of  the  science  instruments.  'Fliis  maneuver  also  leads  wire  boom  interaction  with  a  frequency  of 
about  10  Hz.  The  motion  is  almost  entirely  in  the  spin  plane  with  the  spin  axis  direction  unchanged. 
Since  this  oscillation  is  short  compared  to  natural  spin  rate  changes,  it  does  not  create  particularly  sever 
difficulties  for  attitude  modeling. 


6.4  CANISTER  RELEASES 


The  effect  of  the  release  of  canisters  from  the  spacecraft  depends  on  the  type  and  location.  Single 
canister  releases  generally  give  rise  to  nutation  which  is  much  like  that  following  attitude  adjustments. 
Duel  releases  do  not  in  general  change  the  spin  axis  position.  All  releases  result  in  spin  rate  changes, 
either  positive  or  negative,  which  have  much  the  same  general  effects  as  the  spin  rate  changes  induced 
by  the  thrusters. 


6.5  ORBIT  ADJUSTS 


Orbit  adjustments  to  date  have  been  performed  at  perigee  and  have  been  for  the  purpose  of  raising  the 
altitude  of  apogee.  They  cause  a  decrease  in  the  spin  rate  of  substantial  magnitude  and  the  same  sort 
9f  wire  boom  interactions  as  do  attitude  adjustments,  though  of  smaller  magnitude. 

'fhese  anomalies  are  labeled  in  the  attitude  model  by  codes  which  are  listed  in  Table  2  below.  Of  those 
listed,  only  the  boom  deployment  has  not  been  used  to  date,  'fhis  is  because  the  20  rpm  spin  rate  of  the 
early  portion  of  the  mission  was  not  appropriate  for  attitude  calculation  and  attitude  before  orbit  58  was 
not  routinely  calculated. 


Table  2.  Mission  Event  Codes 


0 

Normal  Segment 

1 

Spin  Up 

2 

Spin  Down 

4 

Canister  Eject 

8 

Attitude  Adjust 

16 

Boom  Deployment 

32 

Orbit  Adjust 

64 

Begin  Eclipse 

128 

End  Eclipse 

256 

Clock  Jump 

512 

Unspecified  Anomaly 
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APPENDIX.  HISTORY  OF  MISSION  EVENTS 


Below  is  a  list  of  the  attitude  events  that  have  been  encountered  through  Orbit  800.  Numbers  in 
parentheses  after  the  canister  releases  identify  the  canister  released. 


Orbit 

Year  Day  Date 

ur 

Event 

058 

1990  230  18  Aug 

24:04 

Attitude  Adjust 

060 

1990  231  19  Aug 

19:43 

Attitude  Adjust 

061 

1990  232  20  Aug 

01:44 

Clock  Jump 

063 

1990  232  20  Aug 

19:59 

Clock  Jump 

063 

1990  233  21  Aug 

01:43 

Attitude  Adjust 

066 

1990  234  22  Aug 

04:51 

Attitude  Adjust 

067 

1990  234  22  Aug 

16:30 

Attitude  Adjust 

071 

1990  236  24  Aug 

08:53 

Clock  Jump 

075 

1990  237  25  Aug 

16:16 

Clock  Jump 

090 

1990  244  01  Sep 

02:09 

Attitude  Adjust 

091 

1990  244  01  Sep 

13:27 

Attitude  Adjust 

107 

1990  251  08  Sep 

03:28 

Attitude  Adjust 

112 

1990  253  10  Sep 

05:44 

Canister  Eject  (20,44) 

122 

1990  257  14  Sep 

08:19 

Canister  Eject  (18,42) 

129 

1990  260  17  Sep 

03:47 

Clock  Jump 

132 

1990  261  18  Sep 

04:54 

Clock  Jump 

138 

1990  263  20  Sep 

21:33 

Attitude  Adjust 

150 

1990  268  25  Sep 

15:28 

Clock  Jump 

154 

1990  270  27  Sep 

08:25 

Clock  Jump 

157 

1990  271  28  Sep 

13:12 

Qock  Jump 

161 

1990  273  30  Sep 

02:54 

Attitude  Adjust 

164 

1990  274  01  Oct 

05:09 

Attitude  Adjust 

165 

1990  274  01  Oct 

22:10 

Clock  Jump 

187 

1990  284  11  Oct 

01:01 

Attitude  Adjust 

190 

1990  285  12  Oct 

01:32 

Attitude  Adjust 

199 

1990  288  15  Oct 

20:29 

Clock  Jump 

199 

1990  288  15  Oct 

21:21 

Clock  Jump 

209 

1990  293  20  Oct 

01:00 

Clock  Jump 

210 

1990  293  20  Oct 

05:59 

Attitude  Adjust 

212 

1990  294  21  Oct 

06:25 

Attitude  Adjust 

232 

1990  302  29  Oct 

06:56 

Attitude  Adjust 

234 

1990  303  30  Oct 

07:57 

Attitude  Adjust 

241 

1990  306  02  Nov 

03:41 

Qock  Jump 

256 

1990  312  08  Nov 

01:44 

Attitude  Adjust 

258 

1990  312  08  Nov 

21:31 

Attitude  Adjust 

274 

1990  319  15  Nov 

14:12 

Clock  Jump 

278 

1990  321  17  Nov 

09:28 

Attitude  Adjust 

281 

1990  322  18  Nov 

07:48 

Attitude  Adjust 

290 

1990  325  21  Nov 

24:52 

Spin  Down 

A-1 


291 

1990 

326  22  Nov 

11:17 

Spin  Down 

302 

1990 

330  26  Nov 

23:03 

Attitude  Adjust 

305 

1990 

332  28  Nov 

07:38 

Attitude  Adjust 

310 

1990 

334  30  Nov 

06:38 

Qock  Jump 

3i7 

1990 

337  03  Dec 

03:10 

Attitude  Adjust 

322 

1990 

339  05  Dec 

06:25 

Clock  Jump 

336 

1990 

344  10  Dec 

23:14 

Attitude  Adjust 

338 

1990 

345  11  Dec 

21:48 

Attitude  Adjust 

342 

1990 

347  13  Dec 

12:35 

Unspecified  (see  1)* 

345 

1990 

348  14  Dec 

16:04 

Unspecified 

360 

1990 

354  20  Dec 

24:07 

Attitude  Adjust 

363 

1990 

356  22  Dec 

04:04 

Attitude  Adjust 

391 

1991 

002  02  Jan 

16:39 

Qock  Jump 

392 

1991 

003  03  Jan 

03:19 

Attitude  Adjust 

395 

1991 

004  04  Jan 

01:43 

Attitude  Adjust 

417 

1991 

013  13  Jan 

01:50 

Canister  Eject  (16) 

417 

1991 

013  13  Jan 

06:37 

Canister  Eject  (17,41) 

422 

1991 

015  15  Jan 

03:42 

Canister  Eject  (40) 

424 

1991 

016  16  Jan 

05:58 

Canister  Eject  (13) 

429 

1991 

018  18  Jan 

04:53 

Canister  Eject  (24,48) 

434 

1991 

020  20  Jan 

05:03 

Canister  Eject  (21,45) 

446 

1991 

025  25  Jan 

01:18 

Attitude  Adjust 

448 

1991 

026  26  Jan 

01:42 

Attitude  Adjust 

451 

1991 

027  27  Jan 

03:05 

Attitude  Adjust 

453 

1991 

028  28  Jan 

02:12 

Attitude  Adjust 

480 

1991 

039  08  Feb 

02:48 

Attitude  Adjust 

482 

1991 

040  09  Feb 

02:53 

Attitude  Adjust 

490 

1991 

043  12  Feb 

03:49 

Canister  Eject  (22,46) 

502 

1991 

048  17  Feb 

03:03 

Canister  Eject  (19,43) 

524 

1991 

057  26  Feb 

06:28 

Attitude  Adjust 

527 

1991 

058  27  Feb 

06:19 

Attitude  Adjust 

556 

1991 

070  11  Mar 

04:31 

Attitude  Adjust 

558 

1991 

071  12  Mar 

02:37 

Attitude  Adjust 

590 

1991 

084  25  Mar 

02:06 

Qock  Jump 

592 

1991 

085  26  Mar 

04:33 

Qock  Jump 

592 

1991 

086  27  Mar 

04:44 

Qock  Jump 

596 

1991 

086  27  Mar 

14:49 

Qock  Jump 

597 

1991 

087  28  Mar 

04:16 

Qock  Jump 

597 

1991 

088  29  Mar 

05:32 

Qock  Jump 

600 

1991 

088  29  Mar 

04:18 

Qock  Jump 

601 

1991 

088  29  Mar 

13:11 

Qock  Jump 

601 

1991 

088  29  Mar 

14:18 

Qock  Jump 

^  This  and  the  following  anomaly  have  been  traced  to  the  failure  of  one  of  the 
batteries.  This  apparently  caused  a  gas  jet  leading  to  a  gradual  spin  rate  decrease.  The 
time  given  is  the  start  of  the  spin  rate  change. 
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788 

1991 

165 

14 

Jun 

789 

1991 

165 

14 

Jun 

789 

1991 

165 

14 

Jun 

789 

1991 

165 

14 

Jun 

790 

1991 

166 

15 

Jun 

790 

1991 

166 

15 

Jun 

791 

1991 

166 

15 

Jun 

793 

1991 

167 

16 

Jun 

795 

1991 

168 

17 

Jun 

796 

1991 

168 

17 

Jun 

796 

1991 

169 

18 

Jun 

797 

1991 

169 

18 

Jun 

799 

1991 

169 

18 

Jun 

799 

1991 

170 

19 

Jun 

15:39  Spin  Up 

18:47  Qock  Jump 

21:59  Qock  Jump 

22:10  Qock  Jump 

05:19  Qock  Jump 

13:28  Qock  Jump 

15:21  Qock  Jump 

11:23  Attitude  Adjust 

12:52  Qock  Jump 

22:16  Attitude  Adjust 

02:12  Qock  Jump 

13:11  Qock  Jump 

24:59  Qock  Jump 

01:16  Attitude  Adjust 
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